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Abstract 

The gap anisotropy in Bi2Sr2CaCu208+5 is revisited in the framework of a d-wave 
scenario in view of the recent angle-resolved photoemission experiment. Based on a 
tight-binding fit to the normal state dispersion, a detail analysis on the effects of the 
inclusion of the next harmonic in the d-wave has been presented. Significant effect 
has been observed in the superconducting Tc. The density of states is linear at the 
nodes with enhanced weight, caused by a marked increase in the low energy excitaions 
which affect the thermodynamics considerably. The slope of the ps — T curve in the 
low temperature regime increases and the specific heat reflects the enhanced entropy 
at low temperatures. The leading edge of the ARPES energy distribution curves have 
been calculated and found to shift towards higher energy. The effect of scattering by 
non-magnetic impurities in this context are also outlined. 

PACS Nos. 74.72-h, 74.20.Fg 
Introduction 

The nature of the superconducting gap anisotropy in the high temperature superconduc- 
tors has attracted a lot of attention in the last few years|l|, 0. A considerable progress has 
been made in the recent past with the resolution and clarity of the angle-resolved photoemis- 
sion spectroscopy (ARPES) having reached a very high level backed by careful and thorough 
analysis of the data]^. A general consensus seems to have emerged that the symmetry of the 
order parameter (OP) in most of the high temperature superconductors is predominantly 
d-wave. In a very recent ARPES experiment on the underdoped Bi2212 system Mesot et. 
al.[^ have reported an observable departure from a simple interpretation in terms of the 
usual d-wave symmetry. Inclusion of the next higher harmonic in the superconducting pair- 
ing function seems to give a better fit to their data. We develop the idea in order to look 
for other observable effects of such a term, and calculate its effects on various physical quan- 
tities in the framework of the usual phenomenological theory[0 in the weak coupling limit 
and suggest further experiments. We have calculated the changes in superconducting tran- 
sition temperature, nodal structure on the fermi surface (FS), the density of states (DOS) 
in the superconducting state and its effect on the specific heat. The energy distribution 
curves (EDO) (that are seen in the ARPES experiments) are also obtained. Fluctuations 
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are included in one loop level and the superfluid density (and hence the penetration depth) 
is obtained as a function of temperature. We observe the differences in slopes of Ps(T) in 
the low temperature regime with increasing level of mixing of the higher harmonic. We also 
calculate the effect of both doping and the mixing of higher order terms on the specific heat. 

Model and Calculations 

In the experiments of Mesot et al, ARPES results have been fitted to a d-wave gap 
function for several Bi2212 systems. The maximum gap values have been adjusted at each 
of the doping levels for the best fit, and the EDCs at different angles and the slope of the 
gap function close to the node have been carefully fitted. Clear deviation from the usual 
dx2_y2 behaviour is observed and attempts were made to fit the data with an admixture 
of about 4-10% of the next harmonic cos(60). A much improved fit was indeed obtained 
with the inclusion of this higher order term. As one underdopes the system, the normal 
state resistivity is known to increase and effective screening becomes weaker. Under such a 
situation it is reasonable to assume that the higher harmonics in the effective interactions 
have to be brought in|p, and highly sensitive experiments would reveal the effects due to 
the growing range of interaction. In the spin fluctuation models too the pairing interaction 
grows sharply in the momentum space with underdoping and long range effects in real space 
become increasingly relevant. 

We undertake to examine the ramifications arising out of the addition of the higher 
harmonic with a model band structure that reproduces the observed FS and the location 
of the van- Hove singularity in Bi2212 accurately [^, 0. The six parameters describing the 
hopping on the BiO plane usedg, 0] to obtain the fit were [0.131, -0.149, 0.041, -0.013, -0.014, 
0.013] (in eV), corresponding to a doping of 6 = 0.17. The van-Hove singularity, due to the 
saddle points k = (vr, 0), (0, vr) is 30meV below the FS, as observed in ARPES. 

Superconductivity occurs through attractive long range interactions of which we keep 
only the near-neighbour part^, |^. Combined with the fact that there is a reasonably strong 
on-site repulsion between the electrons in the cuprates, forcing the pairing function \E'(r = 
0) = 0, such an interaction is known to support d-wave pairing|5|, ^, The vanishing of 
the on-site part of the pairing function implies that the J2q'^{^) = 0, where the q-sum 
runs over the entire Brillouin zone (BZ), thereby contributing equal regions of positive and 
negative sign to the OP. The usual d^i-yi (or cos(20)) OP admits of such constraints. All 
the higher harmonics, cos [(2 + 4n)0] are also admissible under the symmetry restrictions 
of the d-wave. In the underdoped systems, it is likely that higher neighbour interactions 
contribute increasingly, but in the present analysis we kept only the first of such terms. We 
also include only the singlet component of the order parameter in our analysis, as there is 
no reason so far to include the triplet part in the parameter regions that one works with in 
these system]^, ^. 

The superconducting gap equation that is numerically solved is the usual mean-field 
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factorized gap equation 



where \^(k — k') is written in the separable form[^ 

V(k-k')=^?r/(k)r7(k'). 

Here we have chosen the basis function ri(k) to be the Bi representation of the one dimen- 
sional irreducible representations of C^y, ?7(k) = ^{coskx — cosky) (the usual d^^-yi symme- 
try); g measures the strength of the attractive interaction. The triplet channels of pairing 
have not been considered in the foregoingP|. The value of g was chosen (320 meV) such that 
the transition temperature Tc at 5 = 0.17 remains close to 80K (without any admixture of 
cos (60) term). In the spirit of Mesot et al., we introduced the next order term cos (60) in 
r](k) as rjik) = acos(20) + (1 — q;)cos(60), where a measures the relative contributions of 
the two terms. 

The superconducting gap has been calculated for each filling (5) at different levels of 
mixing of the higher harmonic. As representative curves, we show in Fig.l the ones with 0, 
4, 10 and 20% of mixing for 6 = 0.17. The solutions of the gap function for different a show 
typical square root behaviour as T — > Tc. The transition temperature (and A(0)) reduces 
considerably as a deviates from one (the inset to Fig. 1). The gap function with the same 
set of values of a are drawn in the Y-quadrant (Fig. 2) for a demonstration of how they 
change with the mixing of cos(60). Similar shifts have been seen by Mesot et al. as well. 

Calculation of the density of states (DOS) in the superconducting state with and without 
the higher harmonic term is straightforward. For a = 1, one gets the usual d-wave DOS 
with p{E) oc \E\ for E ^ 0. Addition of the higher harmonic term makes the rise steeper as 
there are now more excitations available at lower energy though the low energy behaviour 
of the DOS remains the same (linear). In Fig. 3 the DOS is shown for the normal and the 
broken symmetry states (with only a = 1, as the a = 0.8 curve is indistinguishable in the 
scale of that figure). The typical d-wave V-shaped DOS is obtained with the shoulders at 
±A around E = 0. The inset shows for a comparison, in an enhanced scale of energy, the 
effect of the addition of the higher harmonic term. The pile-up of states at lower energy is 
quite evident. 

In order for an analysis of the ARPES spectra we take the point of view|I^ that the 
spectral function 74(k, cu), convoluted with the fermi function f{uj) and the resolution of the 
detector, gives the EDO. The ARPES intensity (without the detector resolution) is given 
by Iof{uj)A(k,Lj) where, Iq is a prefactor weakly dependent on T and u. It depends on 
the incident photon energy, momentum and the (electron-phonon) matrix element between 
the initial and final states. The momentum resolution of the detectors used in experiments 
considered is about one degree in the BZ and is assumed constant over a circular window 
of 1° radius. With the incident photon energy around 20eV, the energy resolution of the 
detector is taken to be a Gaussian of standard deviation 7meV, a value consistent with the 
present day experimental resolutions[p!0|]. 
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With all these taken into account, the ARPES intensity is given by 

/(k,a;) = /o / / Giuj-u)fiu)Aik',u). 

Jk' Ju 

Here G{uj) is the Gaussian energy resolution function discussed above and the k'-integration 
is within a circular radius of 1''. The spectral function is the usual mean-field one 

A(k, u) = ul6{uj - Ek) + vl6{uj + E^) 

where, E]^ = y^e^ + A^, and the coherence factors = |(1 + ek/Ek),v^ = ^(1 — e^/Ei^). 

Taking the normal state dispersion discussed above, A(k, a;) at five representative an- 
gles (on the FS) = 0, 10, 23, 35, 45 have been calculated and their frequency and momentum- 
resolved behaviour (Fig. 4) obtained. The angles are measured with respect to the line 
(tt, tt) — (tt, 0) as is the practice[0]. We have restricted ourselves to the Y-quadrant of the 
Brillouin zone in order to avoid the complications with the shadow bands that appear in the 
X-quadrant due to the presence of an incommensurate superlattice (along T — Y direction) in 
the BiO planes|^, 0]. We choose a point on the FS (at the above angles), and perform the res- 
olution averaging around that point in momentum and frequency. Extreme care is required 
in evaluating the averages on the FS where the band is most dispersive, i.e., around 45°. 
Panel 1 in Fig. 4 shows frequency averaged A(k, tu) and panel 2, the momentum-averaged 
one. The procedure has been repeated for a = 1.0 and 0.90 to show the effects of the higher 
harmonic on the gap function as one moves along the FS. 

With the momentum and frequency averaged A{k, u) obtained, it is easy to calculate the 
ARPES intensity around different points on the FS. The curves at different angles are shown 
in panel 3 of Fig. 4. The maximum gap is taken to be 30 meV, the same as used by Mesot 
et al, for demonstration. The temperature used is lOK, a typical value used in experiments 
at low temperatures. 

It easily follows from the above discussions that the thermodynamic properties will be 
affected by the addition of the higher harmonics. We calculate the specific heat as a function 
of temperature (Fig. 5) in the superconducting state by taking a derivative of the entropy 
with temperature. Figs. 5 (a), (b) and (c) show the specific heat for different values of a at 
several levels of doping. It is clear from the graphs (a) - (c) that in order to account for the 
extra entropy, the curves move up as higher harmonics are brought in and then the area is 
conserved with reduced transition temperatures. As the doping increases, the curvatures of 
the Fermi surface affect the specific heat as seen in the normal state specific heat curve in 
(c). At Tc the specific heat jumps to its normal state value as is typical of a second order 
transition. 

The calculation of superfluid density ps has been performed using the standard techniques 
of many body theory^, |Tl|. The diamagnetic and paramagnetic contributions to the current 
(and hence the phase stiffness) are calculated in the linear response by first making a Peierls 
substitution tij = Ujexp^ie / ch) J^^ A.dl in the hopping matrix element. The paramagnetic 
contribution at long wavelengths (via excitations above the condensate) in the linear response 
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theory is obtained in terms of the correlation function 



jpara(q) = Um,^,Um^^, J rfr^^(r)e^-- ( [jf™ (q, ^) , jT '^(-Q, 0)] ) A.(q) . 

The correlation function on the RHS is calculated at the one-loop level. Calculation of 
the diamagnetic contribution (from the Meissner condensate) is straightforward: jdia(Q) = 
~]vl!^ ^k,o-(Ck^a-'^k,(7) Ax(q). In the resulting expression, the OP values were taken as their 
mean-field unrenormalized (by the fluctuations) ones, a procedure that is known to work|Q 
except very close to Tc- The resulting ps —T curves are shown in Figs. 6 for a = 1, 0.9 and 
0.8 at different doping levels. Note that the calculations were done in the gauge Ay = Q and 
the gauge invariance is restored if vertex corrections are included. Such an approximation 
entails neglecting the vortex-like fluctuations in the 2D model[]TT[. The above expression 
for Ps is derived for an isotropic order parameter but is expected to work quite well even 
in the anisotropic case at hand as the asymptotic form of the vortex-vortex interaction (at 
high vortex density) is logarithmic and corrections to the expression above due to these 



fluctuations are indeed small. 

Results and Discussion 

The nature of the OP gleaned from Fig. 1 shows a sensitive dependence on a, going down 
as mixing increases. The cos(60) term changes sign four times now in each quadrant and 
such rapid changes average out to a smaller value |P. This is more pronounced if the next 
higher harmonic (cos(lO0)) is introduced [|l^]. As a decreases, the A versus curve becomes 
flatter around the node (Fig. 2) enhancing the quasiparticle excitations above the condensate. 
Such excitations will reflect in a reduced Tc (and A(0)) and pile up of states in the DOS 
at low energies (Fig. 3). This will, of course, affect the thermodynamics considerably. For 
example, though the speciflc heat follows the typical (mean-field) T^ behaviour of a 



d-wave superconductor at low temperature, the increased entropy at lower energy (for a 
deviating from one) manifests itself in the speciflc heat curves (Fig. 5). 

The ARPES line shapes are shown in Fig. 4. The frequency and momentum- averaged 
spectral function is plotted in panel 1 and panel 2 at different angles relative to (tt, tt) 
on the FS. A{\i,uj) (not shown) has the usual delta function behaviour, two sharp peaks 
separated by 2 A at = and closing in as cf) increases and flnally merging at ^ 45°. The 
frequency broadening (panel 1) is a mere consequence of the Gaussian averaging procedure 
on the A(ki?, uj). Strong angle dependence is observed in panel 2 where momentum averaged 
A(^p,uj) is plotted for different 0. The OP does not change within the k-window for small 
angles, while for large angles, the effects are very strong. At large angles the band is highly 
dispersive and large changes in energy occur within the k— window. The k-dependence of the 
OP further enhances this angle dependence since close to the node Ak varies linearly with k 
while at the gap-maximum the variation is a weaker quadratic one. The k-averaged spectral 
weight for a = 0.9 shows perceptibly smaller gap as (p increases. Such a reduction of the gap 
has been observed by Mesot et al. The effect of a is clearly visible in the EDCs in panel 3. 
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The leading edges for the curves corresponding to a = 1 and 0.90 move continuously away 
from the fermi energy as decreases, the curve for a = 1 moving more rapidly. This is the 
scenario depicted in Mesot et al. for seven different angles, where the best fit is obtained 
with 0.89 < a < 0.96 for different samples. The peaks observed in our EDCs are due to the 
mean-field pile-up of states, and not due to electronic correlations [0]. 

The superfiuid density ps is proportional to A^^ (A is the penetration depth) and has 
a power law dependence on temperature jTsU at low temperatures. Figs. 6 (a)-(c) show 
Ps for three representative a (1.0, 0.90 and 0.80) and is seen to fall off faster with the 



inclusion of the higher harmonics. As observed earlier |I5|, |16[, the curves are linear to a 



high degree close to zero temperature. We observe that the slope of ps{T 0) decreases 
monotonically with increasing a (i.e., increasing T^). As temperature rises, excitations 
from the condensate tend to decrease ps at the expense of normal quasiparticles above 
the condensate. The gradual flattening of the A — curve with decreasing a around the 
node makes quasiparticle excitations more accessible at lower temperatures, causing a faster 
descent of ps with temperature. We have also calculated ps at different doping levels and 
find that the slope increases as one underdopes, consistent with the observations of Mesot 
et al. 

It is known that non-magnetic impurities act as pair breaking scatterers for a d-wave 
superconductor [1^ . Repeated scattering even with small momentum transfer on the fermi 
surface between lobes of opposite sign in the Brillouin zone effectively reduces the average 
gap value[]19[ thereby reducing the Tc. Owing to the presence of the cos(60) term, such 
pair-breaking processes are clearly going to be more efficient and will reduce the transition 
temperature as the higher harmonic components increase [|I^]. The pair breaking effect of 
momentum dependent impurity potential and anisotropic gap including the higher harmonics 
have been studied recently |2D|] by solving the Abrikosov-Gorkov equations in the T-matrix 
representation where it has been found that states begin to appear in the gap as the impurity 
potential as well as the higher harmonic component is increased thereby reducing both the 
superconducting gap and the transition temperature. 

In conclusion, we have worked out the effects due to the increasing presence of higher 
harmonics observed recently in a d-wave superconductor on underdoping. Based on the 
standard phenomenological theory for d-wave superconductors!^, we observe that the 
transition temperature, density of states, specific heat, superfiuid density, ARPES intensity 
and the slope of the OP at the node bear clear and detectable signature of this higher 
order term. Although the interaction between quasiparticles have not been included in the 
above, the conclusions drawn remain valid qualitatively on strong physical grounds as has 
been shown in a number of occasions earlier [^, ^. The predictions made here are easily 
verifiable experimentally and will shed light on the nature and strength of the higher order 
term claimed to be present in these superconductors. More experiments are also required to 
fully understand the origin and the physics behind such additional terms. 
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Figure captions 



Fig. 1. Normalized (by A(0) at a = 1) gap versus temperature for four different values of a. 
The inset shows variation of the transition temperature with a. 

Fig. 2. The order parameter as a function of angle (measured with respect to the line (tt, tt) — 
(tt, 0) in the first BZ as in refs.|Q, |^) with increasing mixture of the higher harmonic. 

Fig. 3 The Density of States in the superconducting state (dotted line) with a = 1 and in 
the normal state (solid line) at 5 = 0.17. The inset shows the difference in the DOS 
for a = 1 and a = 0.8 in the superconducting state close to the node at very low 
temperature (the gaps were 15 meV and 11.8 meV as in Fig. 1). 

Fig. 4. The frequency- and momentum-averaged A(k, uj) (panel 1 & 2) and the corresponding 
EDCs (panel 3, see text) at various angles (measured from Y — M direction) on the 
FS (shown outside the panel in the Y-quadrant). Figures (a)-(e) are for angles 45, 35, 
23, 10 and degree. The solid and dotted lines correspond to a = 1.0 and 0.9. 

Fig. 5 The specific heat curves at three different levels of mixing of the higher harmonic ((a), 
(b) and (c) correspond to 6 = 0.10,0.17 and 0.28). The figures clearly reveal the 
enhanced quasiparticle excitations as a deviates from one. 

Fig. 6. The superfluid density is shown against temperature for three different values of 5 = 
0.10, 0.17 and 0.28 (Figs, (a)-(c)). The change in slope at low temperatures is clearly 
visible. 
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